Domain Wall Depinning in Random Media by AC Fields 
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The viscous motion of an interface driven by an ac external field of frequency Uo in a random 
medium is considered here for the first time. The velocity exhibits a smeared depinning transition 
showing a double hysteresis which is absent in the adiabatic case uio ~ > 0. Using scaling arguments 
and an approximate renormalization group calculation we explain the main characteristics of the 
hysteresis loop. In the low frequency limit these can be expressed in terms of the depinning threshold 
and the critical exponents of the adiabatic case. 
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PACS numbers: 75.60.-d, 74.60.Ge 



The driven viscous motion of an interface in a medium 
with random pinning forces is one of the paradigms of 
condensed matter physics |y, |2[ ||. This problem arises, 
e.g., in the domain wall motion of magnetically or struc- 
turally ordered systems with impurities [Q or when an in- 
terface between two immiscible fluids is pushed through 
a porous medium Closely related problems are the 
motion of a vortex line in an impure superconductor Q , 
of a dislocation line in a solid |Q or driven charge den- 
sity waves ||. For a constant external driving force this 
problem has been considered close to the zero tempera- 
ture critical depinning threshold || [u]. [Ll| [l2| and in the 
creep region |jj [l3| . More recently, the results of this ap- 
proach have been used to study hysteresis effects in mag- 
nets subjected to an external force changing adiabatically 
in time @ [n§. It is the aim of this paper to develop 
a description of pinning phenomena in an ac-field in the 
weak pinning limit. As a main result we find that the zero 
temperature depinning transition is smeared and shows a 
pronounced velocity hysteresis. The latter has to be dis- 
tinguished from the hysteresis of the magnetization which 
persists also in the adiabatic case 0, [l5j . Using numer- 
ical simulations, scaling arguments and an approximate 
renormalization group (RG) calculation we derive scaling 
laws for the velocity hysteresis at zero temperature. Also, 
we discuss the influence of thermal fluctuations briefly. 
Despite the fact, that we use the terminology of driven 
interfaces, the results apply correspondingly also to all 
the other systems mentioned above. 

Model and zero frequency critical depinning. — We fo- 
cus on a simple realization of the problem, the motion 
of a D-dimensional interface profile z(x, t) obeying the 
following equation of motion [ fl6[ 

1 dz 

-—= TV 2 z + h cos cj t + g(x,z). (1) 
7 at 

7 and Y denote the mobility and the stiffness constant 
of the interface, respectively, and h(t) = hosmtoot is 
the ac driving force. The random force g(x, z) is as- 



sumed to be Gaussian distributed with (g\ — and 
(g(x, z)g(x\ z')) = S D (x - x')A (z - z'). We further as- 
sume Ao(z) = Ao(— z) to be a monotonically decreasing 
function of z for z > which decays to zero over a finite 
distance I. Under these conditions the relation between 
applied force h{t) and the average velocity (z) = v shows 
in the steady state the inversion symmetry h — > — h, 
v —>■ —v (cf. Fig. [I]). We therefore restrict ourselves 
to the region h > in the further discussion. 

In eq. (|l|) was considered in the adiabatic limit 

luq — > 0. In this case the interface undergoes a second 
order depinning transition at ho = hp where the velocity 
v vanishes as a power law v ~ (ho — hp)@ for ho \ hp, 
(3 < 1. At ho — hp the interface is self-similar with 
a roughness exponent Cj < £ < 1, and the dynamics 
is superdiffusive with a dynamical exponent z, 1 < z < 
2. The critical exponents were calculated up to order 
e = 4- Din|||o| and recently to order e 2 in 
and are related by the scaling laws f3 — v(z — £) and 
v = 1/(2— C) M, where v denotes the correlation length 
exponent: £o ~ \ho — hp\~ v . For ho / hp the divergence 
of £o is related to the increasing size of avalanches. 

In the case of an ac-drive the behavior of the system is 
governed by the two dimensionless quantities ho /hp and 
Cc'o I uj p , where uip = jhp/l. In this paper we mainly focus 
on the most interesting case < u>o <C wp and ho > hp 
since this is the region where universality is expected to 
hold. As illustrated by the numerical solution of eq. (|l|) 
for D = 1 at finite ujo, the sharp depinning transition 
is replaced by a velocity hysteresis, which has clockwise 
rotation: the velocity reaches zero at h(t) = ±h c for 
decreasing and increasing field, respectively (see Fig. [l|). 
For ho 3> hp a second weak hysteresis is found in the 
region h > h c w hp which has anticlockwise rotation. 

Scaling considerations. — First we consider the rele- 
vant length scales of the problem, beginning with lu = 0. 
Comparison of the curvature and the random force term 
on r.h.s. of (|l|) shows, that weak random forces accumu- 
late only on the Larkin scale L P ss ((r0 2 /A (0)) 1/(4 ~ C ' ) 
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FIG. 1: Numerical solution of eq. (y) for a ID interface in the 
dc (solid line) and ac case (dotted line) for ui — wo/wp = 0.3 
and ho > hp ~ 0.27. x is discretized in N = 1000 sites 
and g random in [—0.5,0.5]. The solution is averaged over 
typically 100 disorder (g-) configurations. The arrows show 
the direction of the hysteresis, for \h(t)\ < hp it is clockwise 
and anticlockwise for \h(t)\ > hp. The dashed line shows the 
hysteresis in the case ho ~ hp. 



to a value comparable to the curvature force. On scales 
L < Lp the curvature force density TIL~ 2 is larger than 
the pinning force density, the interface is essentially flat 
and hence there is no pinning. For L > Lp pinning 
force densities exceed the curvature forces, the interface 
becomes rough and adapts to the spatial distribution of 
pinning forces. The largest pinning force density then 
results from L w Lp from which one estimates the de- 
pinning threshold h p w WLp 2 . On scales L ^> Lp 
perturbation theory breaks down. The RG calculation 
performed in g |l(| fy], Q resulted in a scale depen- 
dent mobility and renormalized pinning forces. A finite 
frequency ujq of the driving force acts as an infrared cut- 
off for the propagation of perturbations, resulting from 
the local action of pinning centers on the interface. As 
follows from (|l|) (with 7 — * j(L/Lp) 2 z for L > Lp 
H) these perturbations can propagate up to a length 
scale L^ = L P ( 7 r /uj L 2 p ) 1 / z = L p {u;p/luo) 1/z @: (i) 
If L u < Lp, i.e., loq > top, z has to be replaced by 2. 
During one cycle of the ac drive, perturbations resulting 
from local pinning centers affect the interface configura- 
tion only up to scale L w , such that the resulting curvature 
force is always larger than the pinning force - there is no 
pinning anymore and the velocity hysteresis disappears. 
(ii) In the opposite case L u > Lp, i.e., ujq < u)p, the 
pinning forces can compensate the curvature forces at 
length scales larger than Lp. As a result of the adaption 
of the interface to the disorder pinning forces are renor- 
malized. This renormalization is truncated at L w . In the 
following we will argue, that, contrary to the adiabatic 
limit luq — > 0, there is no depinning transition if ujq > 0. 
Indeed, a necessary condition for the existence of a sharp 
transition in the adiabatic case was the requirement, that 
the fluctuations of the depinning threshold in a corre- 



lated volume Shp « /ip(£p/£o)' D+ ^ 2 are smaller than 
(h-hp), i.e., {D + Qv > 2 g. For uj > the correlated 
volume has a maximal size L u and hence the fluctuations 
Shp are given by 



Shp 
hp 



(D+O/2 



%L) . (2) 



Thus, different parts of the interface see different de- 




transition 



en 



smeared 



irO) adiabatic 
''^T region 



h 



FIG. 2: Schematic frequency-field diagram for the depinning 
in an ac external field (with ho > hp): For < wrj -C up 
the depinning transition is smeared but traces of the luq = 
transition are seen in the frequency dependency of the velocity 
at h — hp. This feature disappears for wo > wp. 

pinning thresholds - the depinning transition is smeared. 
Shp has to be considered as a lower bound for this smear- 
ing. A full understanding of the velocity hysteresis re- 
quires the consideration of the coupling between the dif- 
ferent L w -segments of the interface, which we will do 
further below. When approaching the depinning transi- 
tion from sufficiently large fields, ho 3> hp (ujq <C wp), 
one first observes the critical behavior of the adiabatic 
case as long as £0 <§C L u . The equality £ sa L^ defines 
a field h co signaling a cross-over to an inner critical re- 
gion where singularities are truncated by L^. Note that 
h c0 - hp = /i P (w /wp) 1/( " z) > Shp (cf. Fig. g). It is 
then obvious to make the following scaling Ansatz for the 
mean interface velocity (ho > hp, vp = uipl) 



v (h(t)) as Vp ( — — 

Up 
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hp 
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(3) 



Here the subscript ± refers to the cases of h ^ 0, re- 



spectively, and 4>±[x — > 00] 



(for h 



hp the 
<f>± ap- 



h P > 

classical exponent (3=1 applies). For |x| -C 1, 
proaches a constant c±. Function <j>_ changes sign at 
a critical value h c (u)o) w hp(l — ^(luq/ujp) 1 /^^). Di- 
rect numerical solution of eq. (Q) in D = 1 is in good 
agreement with this prediction as shown in Fig. [| 

Renormalized perturbation theory. — To consider the 
coupling between different L^-segments we treat model 
(|l|) in perturbation theory. After going over to a co- 
moving frame one obtains in lowest non-trivial order in 
g the following equation for the velocity v — io(t), with 
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FIG. 3: Numerical results for v(hp), using eq. (|l|), as a func- 
tion of uj = u)q/lup in D — 1. The dashed line shows the 
prediction of eq. (4) with = 0.17, in good agreement 

with the value 0.19 found in ]l^ ]. 



*b(f) = (z(x,t)) @: 

I„(t) = ft( t ) + / (ft 7 / <^ p («o(t))ff_ p (* (* - «0)> x 
7 Jo 
r 7P 2 t' 



x 7e 



/i(t) + r (t). 



(4) 



Here Jp = / (Cp and 3p( z ) = / rf D ^P x .g(x,z). Re- 
placing the pair correlator of the random forces by 
AoOo(M')) where z (t,t') = J*_ t , dt" v(t") we get 



ro(t) 



dt'(l+ujpt')- D/2 A'o{z{t,t')). (5) 



which results in corrections to the driving force and to 
the mobility which are in general non-local in time. If 
we assume that Ao(z) is an analytic function of its ar- 
gument, it is easy to show that for ujq ^> ujp, ro(t) is 
of the order (wp/cjo) 2 and hence small. In this param- 
eter region the pinning potential merely slows down the 
motion of the wall in agreement with the result of our 
scaling considerations. In the opposite case, loq <C wp, 
perturbation theory breaks down. To treat this frequency 
region it is instructive to consider first the case of a dc 
drive, h{t) = ho, where the velocity is constant and hence 
z(t,t') = vt'. The t'-integral (^) leads to a correction of 
the mobility which diverges as (l/v)( 4 ~ D ^ 2 for v — > 
and D < 4. This divergence could be removed by a RG 
treatment developed in [^|, [l(J . As a result of the elimina- 
tion of the Fourier components z p > with |p| < |p'| < Lp 1 
from eq. (0), 7 and Aq(z) are replaced there by the 
renormalized quantities 

7(p)~ 7 (pL P )- 2 +*, (6) 
A p (z) » K-\Tl/L P )^- D -^A* (z(pL P )</l)) . (7) 

A* (x) exhibits a cusp-like singularity at x = which de- 
velops on scales larger than Lp. In particular, A* fx) ~ 
1 - V^27;M + (e - C)z 2 /6 + 0(\x\ 3 ) for |x| « 1 |, 
and A*e" A * = e^ 1 ^ 2 / 6 for x > 1 Q. In this way 
one generates a renormalized equation of motion which 



serves as starting point for a convergent perturbative ex- 
pansion. The replacements (^|), (^) are valid for mo- 
menta < p < Lp 1 where £ denotes here the cor- 
relation length £q °f the zero frequency depinning tran- 
sition. In the spirit of the RG treatment fluctuations 
on scales larger than £ can be neglected since they are 
uncorrelated. To get the lowest order corrections in the 
convergent expansion one has to replace in eq. (0) the 
bare quantities by the renormalized ones. This leads in 
the limit v -> to r = h P A*'(Q + )/(2- () = -h P which 
is the RG result for the threshold value hp. Replacing in 
addition 7 by 7(£~ 1 ) on the l.h.s. of eq. (||), one obtains 
the correct result for the critical behavior of the velocity: 
v « ^/Lpf-^h - hp) « v P ((h - hp)/h P f. 

In the case of an ac drive the velocity v(t) is peri- 
odic with 2tt/luq. In each cycle of h(t) the velocity goes 
through a region of small values, in which perturba- 
tion theory gives a contribution to 7 -1 proportional to 

A-D 

(l/v(t)) 2 as long as the period is large compared to 
l/v(t). The cutoff t c of the i'-integration is given by the 
approximate relation t^ 1 « ujq + v[t)Jl. These contribu- 
tions are still large at u>o,v(t) — > 0. This breakdown of 
perturbation theory can be overcome by using the RG 
results discussed above on intermediate length scales as 
in the dc case. Such a procedure is justified for momenta 
in the range Lp 1 > |p| ^> £ _1 , where £ is the mini- 
mum of £0 and L u . At these scales the interface is still 
at criticality. Since £0 depends via h(t) on time and eq. 
(JI|) includes retardation effects, a time dependent cutoff 
complicates the problem. Therefore we will restrict our 
consideration to the inner critical region where £ ~ L w , 
i.e., \hp — h| < hp(u!o/ujp) 1 ^ uz . The renormalized effec- 
tive equation of motion follows from (0J with the replace- 
ments (|) and (@) @: 



v(t) 



%- = h{t) + hpujp f dt' f dpp 1+z - c 

w ) Jo JLZ 1 

t 



u, P p*t' A *>( dt"v{t")p(/l), (8) 



with L u = [to p I ujq) 1 I ' z and p = pLp. With that form 
of the equation of motion, we can explain the hystere- 
sis appearing for \h\ < hp (cf. Fig. |l|) more detailed. 
First we consider h < 0: At h c the sign of the velocity 
changes although the driving force is still positive. This 
can be understood as follows: until time t = t c , with 
h(t c ) = h c , the velocity was positive during half a pe- 
riod, hence the argument of the A*' function is positive. 
Therefore the second term of the r.h.s of eq. (||) is neg- 
ative and cancels the positive driving force. To solve eq. 
(^) analytically we consider a parameter region where 
the argument of A* (x) is small compared to unity, i.e., 
A*'(x) w A*'(0+)sgn(/ t t _ t , dt"v(t")). One can show a 
posteriori that this condition is satisfied if ho = OQip). 
With this approximation the momentum integral in (ph 



4 



can be calculated, and we get: 



ill- 



h(t) 



hp 

vz 



S(t,LUp) - " S(t,LU ) 



(9) 



Here S(t,ui) = J Q drr Tg(r) sgnzo(i, t/uj), 6 — 

l/(vz) + 1 and f j(t) = r s (0) - T s (t), where T s {t) = 
dtt 5 ~ 1 e~ t . z (t,r/uj) changes its sign at t = to + 
wk/ujo, n £ Z. The dominating part to S(t,u>) comes 
from r < 0(1). To solve this integral equation for h > 
and < wp, we note, that the sign of zo(t,T /u>p) is al- 
ways positive for r < 1 and hence S(t, up) « vz. For 
t < t c also zo(t,r/u>o) > for the dominating small 
r region of the r-integration in S(t, u>o). This leads 
to 4>—(x) ~ C- + i, c_ « S(t c ,Wo)/vz. By decreas- 
ing t, S{t,u>o) is diminished since regions with negative 
Zo(t, t/u>o) contribute increasingly, which in turn explains 
the second weak hysteresis observed in Fig. |l| |2(| . Next 
we consider the region t > t c , i.e., h < h c , v < 0. By 
increasing t, S(t,u>o) is reduced with respect to S(t c , u>o) 
which leads to a positive curvature of v(h) for h < 0. 
Although that region is beyond the scope of our RG cal- 
culation, since retardation effects require to consider the 
avalanche motion in the region h < h c , it is then tempt- 
ing to conclude \v(h = 0,h < 0)| = 0{uj /io P f/ uz . For 
large negative values of h, v(h) has to reach again the 
result of the adiabatic limit. Together with the inversion 
symmetry this explains the inner hysteresis. 

Thermal fluctuations. — Finally we consider the influ- 
ence of thermal fluctuations on the force - velocity re- 
lation, restricting ourselves to the low temperature re- 
gion T <C Tp = Tl 2 Lp, Tp is a typical pinning en- 
ergy, (i) In the adiabatic limit and for \ho — hp\ <C hp, 
the velocity obeys the scaling relation v(h,T) = (h — 
hp)@ip [(h — hp) 6 /T] where 9 is a new exponent which 
depends on the shape of the potential at the scale Lp 
pl[ . For uj > one can extend the scaling relation (^) 



to a second scaling field h h ^ 1 
particular for h w hp 



v{h P ,T) w v P 



wp 



'T P < 



T 



and one finds in 



WO 



(10) 



with <f>± [x —* oo] ~ x' 3 and [x — > 0] ~ c± . The 
thermal smearing of the zero frequency depinning tran- 
sition is still seen at finite loq as long as u>o < ljtQip) ~ 

J . On the other hand for small fields, h <C 
hp, the domain wall shows creep behavior v(h,T) « 
wpe T v'w , [i = 2 ^^ D r 2 ; where £ denoted the equi- 
librium roughness exponent Q . fiij It was shown in ]15[ , 

that the creep law is valid also at finite frequencies as 

_Tp_/h P y 

lon g as ujq <C uJrjti) ~ tjpe r V fe / , fa -c Zip- For 



- velocity relation is that of the adiabatic case at finite 
temperature. 

To conclude, we have shown, that the sharp depinning 
transition of an interface driven by a dc field is smeared 
showing a pronounced velocity hysteresis, when the ex- 
ternal drive is oscillating. The size of the hysteresis is 
described by the power laws eqs. (||) and (|l^) which 
are supported by an approximate renormalization group 
analysis and a numerical simulation. The case ho < hp 
will be the subject of forthcoming studies. 
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